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Gravitational-wave interferometers are expected to monitor the last three minutes of inspiral 
and final coalescence of neutron star and black hole binaries at distances approaching cosmological, 
where the event rate may be many per year. Because the binary's accumulated orbital phase can be 
measured to a fractional accuracy -C 10~ 3 and relativistic effects are large, the waveforms will be 
far more complex, carry more information, and be far harder to model theoretically than has been 
expected. Theorists must begin now to lay a foundation for extracting the waves' information. 

PACS numbers: 04.30,+x, 04.80,+z, 97.60.Jd, 97.60.Lf 



A network of gravitational-wave interferometers (the 
American LIGO §}, the French/Italian VIRGO § and 
possibly others) is expected to be operating by the end 
of the 1990s. The most promising waves for this network 
are from the inspiral and coalescence of neutron star and 
black hole binaries ||Q , with an estimated event rate of 
- (3/year)(distance/200 Mpc) 3 g. 

A binary's inspiral and coalescence will produce two 
gravitational waveforms, one for each polarization. By 
cross-correlating the outputs of 3 or more interferome- 
ters, both waveforms can be monitored and the source's 
direction can be determined to within ~ 1 degree p]|| . 

We shall divide each waveform into two parts: the in- 
spiral waveform, emitted before tidal distortions become 
noticeable (3 or fewer orbital cycles before complete dis- 
ruption or merger [Q), and the coalescence waveform, 
emitted during distortion, disruption, and/or merger. 

As the binary, driven by gravitational radiation reac- 
tion, spirals together, its inspiral waveform sweeps up- 
ward in frequency /. The interferometers will observe 
the last several thousand cycles of inspiral (from / ~ 10 
Hz to ~ 1000 Hz), followed by coalescence. 

Theoretical calculations of the waveforms are generally 
made using the Post-Newtonian ("PN") approximation 
to general relativity. Previous calculations have focused 
on the Newtonian-order waveforms and on PN 

modulations of their amplitude and frequency ||] . 

Two of us (Tl| recently realized that the PN modula- 
tions are far less important than PN contributions to the 
secular growth of the waves' phase $ = 2ir J fdt (which 
arise, largely, from PN corrections to radiation reaction.) 
The binary's parameters are determined by integrating 
the observed (noisy) signal against theoretical templates 
, and if the signal and template lose phase with each 
other by one cycle (out of thousands) as the signal sweeps 
through the interferometers' band, their overlap integral 
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will be strongly reduced. Correspondingly, one can infer 
each of the system's parameters Xi to an accuracy that is 
roughly the change AAi which alters by unity the number 
of cycles Af cy c spent in the interferometers' band. 

The only parameters Xi that can influence the inspiral 
template significantly are (i) the initial orbital elements, 

(ii) the bodies' masses and spin angular momenta, and 

(iii) if the spins are very large, the bodies' spin-induced 
quadrupole moments. (The quadrupole moments pro- 
duce orbital phase shifts a little larger than unity at 
/ ~ 10 Hz, but negligible phase shifts later j7|. Below we 
shall ignore this tiny effect.) 

If (as almost always is the case) the binary's orbit has 
been circularized by radiation reaction, then the num- 
ber of cycles spent in a logarithmic interval of frequency, 
dMcyc/dlnf = (l/27r)(d$/tfln/), is: 

tWcyc 5 1 f /743 11 (1 N 

din/ 967r^M 2 / 3 (7r/) 5 / 3 { \336 4 M / 

- [4tt + S.O.ja: 1 - 8 + [S.S.ja: 2 + 0(x 2 5 )} . (0.1) 

Here M is the binary's total mass, /x its reduced mass, 
and x = (7rM/) 2 / 3 ~ M/D the PN expansion parameter 
(with D the bodies' separation and c = G = 1). The PN 
correction [0(x) term] is from In the P 15 N correc- 
tion, the 4ir is created by the waves' interaction with the 
binary's monopolar gravitational field as they propagate 
from the near zone to the radiation zone [|l3| , and the 
"S.O." denotes contributions due to spin-orbit coupling 
fl4j| . In the P 2 N correction the "S.S." includes spin-spin 
coupling effects Q plus an expression quadratic in \xjM. 
Since the leading-order, Newtonian contribution to 



Eq. (3.1) gives several thousand cycles of oscillation in 
the interferometers' band, to be measureable the higher- 
order corrections need only be as large as a part in several 
thousand, even when the amplitude signal to noise ratio 
has only its typical value of 8. It is this that gives the 
waves' phasing its high potential accuracy. 

The determination of the binary's parameters (masses, 
spins, orbital elements) is ma de p ossible by the various 
frequency dependences in Eq. ( |0.l| ). Analytic []l5 16 and 



Monte Carlo |17y calculations show that (i) the "chirp 
mass" M c = (]?> 5 M 2 / 5 [which governs the Newtonian 
part of ( p.l| )] will typically be measured to 0.1 per cent, 
and (ii) if we somehow knew that the spins were small, 
then the reduced mass fi would be measured to ~ 1 per 
cent for NS/NS and NS/BH binaries, and ~ 3 per cent for 
BH/BH binaries. (Here and below NS means a ~ 1.4M Q 
neutron star and BH means a ~ 10M Q black hole.) Un- 
fortunately, t he fr equency (x) dependences of the various 



terms in Eq. (0.1) are not sufficiently different to give a 



clean separation between fi and the spins. Preliminary 



estimates |17| |, in which the S.O. term in Eq. (3.1) was 
taken into account but not the S.S. term, suggest that 
the spin//x correlation will worsen the typical accuracy 
of fj, by a factor ~ 20, to ~ 20 per cent for NS/NS, 
- 25 per cent for NS/BH, and - 50 per cent for BH/BH. 
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These worsened accuracies might be improved signifi- 
cantly, however, by the waveform modulations |{J — most 
especially, modulations caused by spin-induced preces- 
sion of the orbit [H . Much additional theoretical work 
is needed to pin down the measurement accuracies for \x 
and the spins. 

The highest accuracy parameter information will come 
from low frequencies where most of the waves' cycles 
are located |l~5|-|l7|]. (It thus is important that the in- 
terferometers achieve good low- frequency performance). 
When the waves have swept to higher frequencies, the 
binary's parameters may be fairly well known, and the 
subsequent, highly relativistic waveforms might be used 
to map out the binary's innermost spacetime geometry 
near the orbital plane. It would be useful to develop 
a parametrized, "PPN-like" multi-theory formalism for 
translating the measured waveforms into spacetime maps 
and for comparing those maps with the predictions of 
various gravitation theories. 

To make optimal use of the interferometers' data 
will require general-relativity-based waveform templates 
whose phasing is correct to within a half cycle or so dur- 
ing the entire frequency sweep from ~ 10 Hz to ~ 1000 
Hz. By examining an idealized limit (fi <C M and vanish- 
ing spins), several of us [ p~3|Jl~S| ] have discovered that, to 
compute the templates with the desired accuracy via PN 
methods will be very difficult. We calculated the waves 
from such a binary to high accuracy using the Teukolsky 
p8[ black-hole perturbation formalism and then fit a PN 
expansion to the results, to obtain fll9|| 



dKyc = PdE/df = J^M//j 
dhxf dE/dt 96tt a; 2 - 5 X 

1 _ 3 81 2 _ 675 3 



1 - ^§-x + Airx 1 - 5 - 4.9a; 2 - 38a; 2 - 5 + 170a; 3 



(0.2) 



Here dE/df is half the ratio of the changes of orbital en- 
ergy dE and orbital frequency d(f/2) as the orbit shrinks, 
and dE / dt is the power carried off by gravitational waves. 
The coefficients in the denominator after the 4-7T were ob- 
tained from the fit; all other coefficients are known ana- 
lytically. The coefficient 4.9 has accuracy ~ 2% ; the 38, 
- 10% ; and the 170, ~ 25%. 

We can use Eq. ([l^) to estimate how the phase ac- 
curacy depends on the order to which the PN expansion 
is carried. If the template were com pute d through P 2 N 
order [all a; 2 ' 5 and a; 3 terms in Eq. ( |0.2| ) omitted], then 
the phase error typically would reach a half cycle when 
the frequency had swept from 10 Hz to only ~ 15 Hz, 
and the total phase error from 10 Hz to 1000 Hz would 
typically be ~ 6 cycles. For interferometers whose peak 
sensitivity is at 70 Hz [Q, one could use such templates 
for a maximum range of ~ 50 Hz to ~ 90 Hz without 
accumulating phase errors larger than a half cycle. If 
the template were computed through P 2 5 N order, there 
would be no substantial improvement. It is not at all 
clear how far beyond P 2 5 N the template must be carried 
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to keep the total phase error below a half cycle over the 
entire range from ~ 10 Hz to ~ 1000 Hz. 

This slow convergence of the PN expansion might be 
improved by cleverness in the way one expresses the ex- 
pansion (e.g., by the use of Pade approximates). How- 
ever, even with great cleverness and fortitude, PN tem- 
plates might never cover the entire inspiral range, from 
~ 10 Hz to ~ 1000 Hz. Evidently, new techniques are 
needed for computing templates to higher accuracy. Two 
such techniques look promising: A "post-Teukolsky ex- 
pansion," and a "weak-reaction expansion." 

The post- Teukolsky expansion would expand in powers 
of fi/M and thus would be useable only for a light body 
orbiting a much heavier black hole. The expansion's first 
step would be the unperturbed hole's Kerr metric; the 
second, the Teukolsky formalism |L8| for the light body's 
first-order perturbations; and the third, incorporation of 
radiation reaction pofl . Each of these first three steps is 
now in hand, though studies of the consequences of the 
radiation reaction are only beginning That the ex- 
pansion must be carried well beyond these three steps is 
evident from the magnitude of the (i/M term in Eq. ( p.l[ ), 
and from recent studies |^] of the influence of a finite 
n/M on the last stable circular orbit. 

The weak-reaction expansion would be a variant of 
numerical relativity in which one expands in powers of 
1/Q = (dE j ' d£)(jt j 'E) _1 (a measure of the effects of ra- 
diation reaction during one orbit). Because 1/Q ~ x 2 ' 5 
(with x the PN parameter) is always 1, this expan- 
sion might produce adequate inspiral templates for all 
binaries. The expansion's first step might be a numerical 
solution of the Einstein equations for the binary's met- 
ric, with standing-wave boundary conditions at infinity 
and at the BH horizon(s) (if any) In co-rotating 

coordinates, the equations would be elliptic and might 
be solvable to high accuracy by relaxation techniques. 
The second step might be to switch from standing-wave 
to outgoing- and downgoing-wave boundary conditions, 
and evaluate the resulting linearized metric perturba- 
tions, and with them the leading contributions to the 
waveforms. The third step might be the leading effects 
of radiation reaction, including the orbital inspiral. To 
obtain templates with the desired phase accuracy, one 
would have to carry the expansion beyond this third step. 

Although an optimally-precise measurement of the bi- 
nary's parameters requires templates faithful to general 
relativity's predictions, a near optimal search for binary- 
inspiral waves can use a family of templates that merely 
span the range of expected waveform behaviors, without 
being closely related to the true general relativistic wave- 
forms [^4| . Theorists need to develop a set of such search 
templates and optimize their efficiency for wave searches. 

Turn, now, from a binary's inspiral waveforms to its 
coalescence waveforms. By the beginning of coalescence, 
the binary's parameters (masses, spins, geometry) will 
be known with fair accuracy, and from its masses the 
nature of its bodies (BH vs NS) should be fairly clear. 
The coalescence waveforms can then be used to probe 
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issues in the physics of gravity and atomic nuclei: 

If the bodies are BH's of comparable mass (and espe- 
cially if their spins are large and not perpendicular to the 
orbital plane) , then the coalescence should produce large- 
amplitude, highly nonlinear vibrations of spacetime cur- 
vature that may reveal aspects of gravity we have never 
seen. Although supercomputer simulations of such coa- 
lescences are being attempted, they are so difficult, es- 
pecially for spinning holes, that they may well not give 
definitive results without observational guidance. Con- 
versely, the observed waveforms will be very hard to in- 
terpret without the guidance of simulations. 

The interferometer measurements discussed above may 
well be achieved in the early years of LIGO/ VIRGO. 
Most of the measurements discussed below are more dif- 
ficult, and may require more mature interferometers. 

For a NS/BH binary with the BH spinning moderately 
fast, the NS should disrupt tidally before plunging into 
the BH. The NS disruption should be quick, as should the 
final coalescence of a NS/NS binary: within about one 
orbit, the NS/BH and NS/NS binaries may get smeared 
into roughly axially symmetric configurations, thereby 
shutting off their waves [^5|,0. The rapid wave shut- 
off and contamination by laser shot noise may prevent 
the coalescence waveforms from revealing more than two 
numbers: the maximum frequency /max reached, and the 
total wave energy £gw in the last ~ 0.02 seconds (when 
the energy emission is maximal). Either number, how- 
ever, would be valuable. It could tell us the NS radius 
-Rns ( or a combination of the two NS radii) [^7| . Since the 
NS masses will already be known from the inspiral wave- 
forms, such measurements on a number of NS's could re- 
veal the NS radius-mass relation i?NS (-^NS ) , from which 
one can infer the equation of state of matter at densities 
from ~ 1 to 10 times nuclear p6[ . 

In preparation for the interferometers' measurements 
of /max and/or EqWi theorists need to model tidal dis- 
ruption in NS/BH binaries and coalescence of NS/NS 
binaries, to determine the dependence of /max and -Eqw 
on i?NS and the binary's other parameters. (Such model- 
ing is also driven by the possibility that these events are 
the sources of observed gamma-ray bursts pi|.) Sim- 
ple arguments |7y27|| suggest that / max will be roughly 
(1/27t)(M ns /-R|s) T ^ 2 , so for a 1.4M NS, a 15 km ra- 
dius will lead to / max ~ 1000 Hz, while a 10 km radius 
will lead to / max ~ 2000 Hz. Similarly, £qw will be pro- 
portional to l/i?NS so a 10 kni NS will emit ~ 50 per 
cent more energy just before disruption than will a 15 
km NS. 

To measure the shutoff frequency / max ~ 1000 Hz in 
the midst of the interferometers' shot noise will probably 
require specially configured ("dual recycled" |29(| ) inter- 
ferometers that operate over adjustable, narrow bands 
A/ around adjustable central frequencies f . Such inter- 
ferometers, operated in concert with broad-band ones, 
could give simple "yes/no" answers as to whether the 
waves swept through their frequency bands. By col- 
lecting such data on a number of binaries, with various 
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choices for the frequencies / , one might zero in on / max 
for various types of binaries, and thence on the NS radius- 
mass relation. The dual-recycled sensitivities required for 
such a program, however, will be difficult to achieve [B7) . 

Even more difficult, but not totally hopeless, will be 
the measurement of £qw via the waves' "Christodoulou 
memory" |3(],^7]] (a slowly changing mean of the wave- 
forms, produced by the 1/r gravitational field of the emit- 
ted gravitons). When the waves' mean is monitored by 
an optimally filtered, broad-band interferometer, it gets 
averaged over a time ~ 0.02 seconds, so the maximum 
level it reaches is proportional to -Egw- Unfortunately, 
to measure the memory's level and thence Bqw with 
even 50 per cent accuracy for binaries at 200 Mpc dis- 
tance is likely to require interferometers modestly better 
than LIGO's so-called "advanced detectors" jjj. 

Turn, next, from a binary's coalescence to the use of its 
waveforms as cosmological probes. Schutz Q has pointed 
out that a binary's distance from Earth (more precisely, 
its "luminosity distance" r£) can be inferred from its in- 
spiral waveforms, and that such gravitationally measured 
distances, when combined with electromagnetically mea- 
sured redshifts z, might determine the Hubble constant 
H a to better precision than heretofore. However, this 
will require identifying electromagnetically the clusters 
of galaxies in which the binaries lie, at least in some sta- 
tistical sense — a difficult task M . 

Remarkably, for NS/BH binaries at z ^ 1, it may be 
possible to determine z as well as tl from the gravita- 
tional waveforms, and then from r^(z) one might infer 
the Universe's Hubble constant H Q , mean density p, and 
cosmological constant A (A variant of the method 

could also work for NS/NS binaries, but NS/BH bina- 
ries emit more strongly and thus can be seen to greater 
distances, making them more promising than NS/NS.) 

The key to measuring z is the expectation (based on 
radio pulsar observations) that NS masses will cluster 
around 1.4M . (If this fails, but there is a sharp cut- 
off in NS masses at some upper limit, then that limiting 
mass can be used as the key to the measurements.) For a 
binary at z > 1 , the waveforms will be measured as func- 
tions of redshifted time, and correspondingly they will 
reveal redshifted masses (1 + z)M^s and (1 + z)Mbu- 
From the measured luminosity distance and redshifted 
masses, it should be fairly clear whether the binary is 
NS/BH or not; and if so, then one can infer a "candidate 
redshift" for the binary, 1 + z can d = (1 + z)M^s/lAM Q . 
If the NS masses cluster around 1.4M©, then a large sam- 
ple of measured distances and candidate redshifts should 
cluster around the true distance-redshift relation. 

This gravitational method of determining cosmological 
parameters might be much more immune to evolution- 
ary and ill-understood systematic effects than are elec- 
tromagnetic methods, because gravitational waves are 
immune to absorption and scattering, and neutron star 
masses might not be sensitive to the Universe's evolution. 

A detailed analysis by one of us |3l| suggests that 
this gravitational approach to cosmology might begin to 



G 



bring useful information when broad-band interferome- 
ters reach the LIGO "advanced detector" sensitivities 
: with one year of observational data, the one-sigma 
accuracies might be SH a ~ 0.01iJ o , 5p ~ 0.1p cl -it, and 
SA ~ 0.3i?o, where p cr it is the critical density to close 
the Universe. 

We conclude by noting that some of the issues dis- 
cussed in this Letter have implications for a possible 
future space-based interferometer called LAGOS, which 
would operate in the band 0.0001-0.03 Hz ||]. We shall 
discuss those implications elsewhere. 
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